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With a variational three-body calculation we study the role of the interplay between the on- 
site Coulomb, Hund's rule, and superexchange interactions on the spinwave excitation spectrum of 
itinerant ferromagnets. We show that correlations between a Fermi sea electron-hole pair and a 
magnon result in a very pronounced zone boundary softening and strong deviations from the Heisen- 
berg spinwave dispersion. We show that this spin dynamics depends sensitively on the Coulomb and 
exchange interactions and discuss its possible relevance to recent experiments in the manganites. 
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The interaction between itinerant carrier spins and lo- 
calized magnetic moments leads to ferromagnetic order in 
a wide variety of systems [l[ . Examples include the man- 
ganese oxides (manganites) Ri-^Aa-MnO^ ( R=La, Pr, 
Nd, Sm, • • • and A= Ca, Ba, Sr, Pb, ■ • • ) 0] and the III- 
Mn-V ferromagnetic semiconductors 0]. Such systems 
are of great current interest due to their novel potential 
applications. For example, the manganites display colos- 
sal magnetoresistance |2j, while ferromagnetic semicon- 
ductors raise the possibility of multifunctional quantum 
devices that combine information processing and storage 
on a single chip with low power consumption [J] . In such 
materials, the magnetic and transport properties are in- 
timately related and can be controlled by varying the 
itinerant carrier concentration and dimensionality. 

In the manganites, n=\-x itinerant electrons per Mn 
atom partially fill a d-band with e g symmetry. Their 
concentration, n, is controlled by the hole doping, x. 
The d-band kinetic energy K is determined by the hop- 
ping energy between the neighboring lattice sites, t ~ 
0.2 — 0.5eV. The itinerant electron spins interact strongly 
with localized spin-S magnetic moments (Hund's rule 
coupling H exch with strength J ~ 2eV> t). S=3/2 
comes from the three electrons in the tightly bound t2 g 
orbitals. This ferromagnetic interaction competes with 
the direct antiferromagnetic interactions (Haf) between 
neighboring local spins, Jaf ~ O.Olt. The largest energy 
scale in the manganites is given by the on-site Coulomb 
repulsion between the itinerant electrons, U ~ 3.5 — 8eV 
(Hu). This Coulomb interaction is generally difficult to 
treat and its effects have received less attention. Here we 
focus on the role of U on the spin dynamics in the con- 
centration range 0.5 < n < 0.8 where metallic behavior is 
observed in both 3D and quasi-2D (layered) manganites. 

The ferromagnetic order in the manganites can be in- 
terpreted to first approximation by invoking the double 
exchange mechanism and the J — > oo limit of the min- 
imal Hamiltonian K + H exc h @, HI- An itinerant car- 
rier is allowed to hop on a lattice site only if its spin is 
parallel to the local spin on that site. The kinetic en- 



ergy is thus reduced when all spins are parallel. This 
favors the ferromagnetic state \F), which describes local 
spins with S z — S on all lattice sites and a Fermi sea of 
spin-| electrons. The above spins are often treated as 
classical, justified for S — ► oo [2j. In this limit, the sys- 
tem can be described by a nearest neighbor Heisenberg 
model with ferromagnetic interaction. Quantum effects 
are often treated perturbatively in 1/S @,13|- To 0(1/5"), 
one thus obtains noninteracting Random Phase Approx- 
imation (RPA) magnons, whose dispersion in the strong 
coupling limit coincides with that of the nearest neigh- 
bor Heisenberg ferromagnet [t|- Such a dispersion was 
observed experimentally for concentrations n > 0.7 [8j. 

However, strong deviations from the short range 
Heisenberg model spinwave dispersion were observed for 
n < 0.7 in both 3D 0, E(| El[ and quasi-2D mangan- 
ites 



12j. Most striking is the strong 



spinwave softening 
This indicates a 



tie strong s 

close to the zone boundary 0, E3, EH 
new spin dynamics in the metallic ferromagnetic phase 
whose physical origin is still unclear [ll[. The proposed 
mechanisms involve orbital degrees of freedom, magnon- 
phonon interactions, disorder, bandstructure effects, and 
the Hubbard repulsion @, 0, M, EI El El • The zone 
boundary softening can be fitted phenomenologically by 
adding lon g ra nge interactions to the Heisenberg Hamil- 
tonian 



MM 



111 found that the above soft- 
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ening increases with x=l-n, while the dispersion for low 
momenta only changes weakly. They argued that none 
of the existing theories can explain these experimental 
trends [111 ]. 

In this paper we study the concentration dependence of 
the spinwave dispersion predicted by the model Hamil- 
tonian H = K + H exch + Hu + H A f @, @, E3 with a 
single e g orbital per lattice site. We treat exactly the 
long-range magnon-Fermi sea pair three-body correla- 
tions induced by the interplay between Hjj and H exc h 
with a variational wavefunction. We show that such cor- 
relations lead to strong deviations from the RPA and 
Heisenberg spinwave dispersions. These deviations, as 
well as the stability of the ferromagnetic order, depend 
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FIG. 1: (Color online) Spinwave dispersion along different di- 
rections ( n=0.6, J =7t, Jaf =0.012t). (a) U=25t (b) U=10t 
(c) Direction r — X: U=10^45t in increments of 5t. (d) 
Deviation from the RPA: 1 — lo/lo rpa at X -point. 



sensitively on Hu- Our approach interpolates between 
the strong/weak coupling and n=0/n.=l limits with the 
same formalism and can therefore address the interme- 
diate interactions and n relevant to the manganites. At 
the same time, it recovers the 1/S expansion Q and ex- 
act numerical results 15, HH as special cases. We find 
that magnon-Fermi sea pair correlations due to U result 
in a pronounced zone boundary spinwave softening that 
increases with x (similar to the experiment [ill ]) m a wav 
that depends on U and J. Our variational calculation sets 
a lower bound on the magnitude of this softening. 

Method — We use the variational wavefunction |Q) = 
Mq\F), where the operator Mq conserves the total mo- 
mentum Q and lowers the z-component of the total spin 
by 1. This spin reversal can be achieved either by lower- 
ing the localized spin z-component, via the collective spin 
operator S~ [3] , or by coherently promoting an electron 



from the spin-f to the spin- J. band; it may also be accom- 
panied by the scattering of Fermi sea pairs. Neglecting 
multipair excitations, the most general Mq is 14 1 



(i) 



where creates a spin-cr, momentum-k electron, v, /i 
(a) label states inside (outside) the Fermi sea. The first 
two terms create a magnon of momentum Q. The last 
two terms describe magnon scattering, Q — > Q + /i — a, 
accompanied by electron scattering across the Fermi sur- 
face, fi — > a (Fermi sea pair shakeup). By setting 
't = $ = we recover the RPA results [lij]. However, 
here the variational parameters Xfi, 4^ and are 



not restricted in any way; unlike in Ref. 17[ , we do not 
assume any particular form or momentum dependence. 
By solving the full variational equations numerically for 
fairly large N xN lattices ( N ~ 20-30), we put an upper 
bound on the spinwave excitation energies wq (with re- 
spect to \F)) that converges with N and thus reflects the 
thermodynamic limit. We can therefore conclude that (i) 
the exact dispersion is at least as soft as our results, (ii) 
wq < means that \F) is not the ground state. 

The wavefunction Eq.{T]) offers several advantages. It 
gives exact results in the two concentration limits n — > 
(one electron) and n — 1 (half-filling). In the special 
cases H{j — Haf — and H e ^, = Haf — it agrees 
very well with exact results [14j, [l(| . Our results also be- 
come exact in the atomic limit t — > [ljj, |l8| and should 
therefore treat local correlations well. While the latter 
dominate in the strong coupling limit, long range corre- 
lations become important as J/ 1 and U/t decrease [3] ■ 
The experiment (ToL fllT | points out the importance of 
long-range interactions. Eq.JT]) treats exactly all correla- 
tions between a single Fermi sea pair and a magnon. The 
only restriction of Eq. ([1]) is that it neglects contributions 
from two or more Fermi sea pairs, which are however 
suppressed for large S [3] and in ID 1(| 1^ . 

Results — FigQ] shows the calculated three-body 
spinwave dispersion for U—25t (FigQJa)) and U=10t 
(FigQTb)). It compares this to the RPA (* = $ = 0) 
and the results of Ref. Q , which we recover by expanding 
the RPA to 0(1/5) and 0[1/(JS + nU)]. Fig^a) also 
compares to the Heisenberg dispersion obtained by tak- 
ing the limit J — > oo, U — of the RPA (rather than by 
fitting). The latter deviates strongly from our interme- 
diate coupling results. While the RPA agrees well with 
Ref. Q , the Fermi sea pair-magnon correlations lead to a 
very strong softening (deviations ^100% from the RPA). 

The on-site Coulomb repulsion U increases the spin- 
wave energies and therefore the stability of the ferro- 
magnetic state \F). Fig. [TJc) demonstrates this hard- 
ening along r — X ((0,0)— >(7r,0)) as U increases in steps 
AU = 5t. While initially the energies increase strongly 
with U, their relative change decreases with increasing 
U. Nevertheless, full convergence to the U — > oo result 
(dashed curve in FigOTc)) om Y occurs for very large U. 

Ref. [l|| treated the effects of strong U by mapping the 
problem to a Hamiltonian with U = [3| and renormal- 
ized hopping t(n). The magnon excitations were then de- 
scribed within the RPA. Due to the increase in the effec- 
tive J/t(n), U resulted in higher spinwave energies. Here 
we show that carrier-magnon correlations beyond the 
RPA, induced by U, lead to a pronounced zone bound- 
ary softening as compared to the RPA. This can be seen 
in Fig. njcl), which shows the percentage deviation from 
the RPA at the X-point as function of U (maximum is 
100%). While the deviations from the RPA decrease with 
increasing U, they remain quite large for the typical U. 

We now focus on the dependence of the X-point en- 
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FIG. 2: (Color online) J 4 (n)/Ji(n) for Jaf = 0.012f, U = 
25i extracted by fitting our results to the fst+4th nearest 
neighbor Heisenberg model. The same behavior is exhibited 
by the spinwave softening compared to the Heisenberg model. 
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FIG. 3: (Color online) (a) Comparison of different approxi- 
mations for D(n) (J=7t, f/=25t), (b) the role of U on D(n), 
(c) Contribution of magnon-pair correlations for different mo- 
menta, (d) Origin of magnon softening. JAF~0-012t. 



ergy on n. Refs.0, [13, II | found that the deviation, 
at this zone boundary, of the nearest-neighbor Heisen- 
berg model dispersion that fits the experiment at small 
Q increases with x=l-rt. The experimental dispersion 
along all directions in the Brillouin zone was fitted 
by a Heisenberg model with both 4th-nearest-neighbor 
( J4) and next-nearest-neighbor ( Ji) exchange couplings; 
2nd- and 3rd-nearest-neighbor interactions were negligi- 
ble [ll|. The ratio J4/J1 oc x becomes strong for n < 0.7 



11] 



Our numerical results can also be fitted very well to 
the J\-Ja Heisenberg model. Fig[2] shows the behavior 
of J&{n)/ Ji(n) (and thus the spinwave softening) for dif- 
ferent J. The crucial role of the pair-magnon correlations 
is clear by comparing to the RPA. The RPA gives small 
J4/J1 (in the strong coupling limit it coincides with the 
nearest-neighbor Heisenberg dispersion |l4(). However, 
the pair-magnon correlations greatly enhance J4/J1 (and 
the softening), typically by a factor 3-4 or higher in Fig. 

J4/J1 increases rapidly with x—l-n until it reaches 
its maximum. For large J/t, J4/ J\ increases more slowly 
with x. This increase is sharp for smaller J, as the ferro- 
magnetic state becomes less stable (compare Figs. Eta) 
and^d)). On the other hand, J4/J1 is small for n >0.7. 

Next we turn to the spinwave dispersion for small Q. 
Its behavior is characterized by the stiffness D(n), ob- 
tained by fitting the small-Q dispersion to the form DQ 2 . 
FigE{a) compares our results to Ref.0], Eq.([T]), and the 
RPA. The pair-magnon correlations decrease D(n) by as 
much as ^100% as compared to Ref.Q and by as much as 
~50% from the RPA. FigOJa) demonstrates a plateau as 
function of n, where D(n) remains fairly constant within 
a wide range of n relevant to the manganites. The pair- 
magnon correlations decrease the dependence of D on 



n for such concentrations (compare the three curves in 
FigGUa)). As shown in FigJ3£b), U increases the stiff- 
ness. Overall, Figs El and [3] are consistent with the main 
experimental trends 13, HI- However, in Ref. 11 1 D(n) 



was found to be fairly constant over a wider range of n. 
FigEJand FigOshow that the pair-magnon correlations 
suppress the dependence of D on n while enhancing J4. 
We speculate that the differences from the experiment 
may be due to the bandstructure effects neglected here. 

We now turn to the origin of the zone boundary soften- 
ing and show that it is dominated by strong correlations 
due to U. We set Jaf — 0. Similar to Ref.pjj], the spin- 
wave dispersion wq is determined by the amplitude Xj*, 
Eq.([T]). describing the coherent spinf — >spin| electron ex- 
citation (oc X^^m) an( l by tue amplitude 'J?, describing 
magnon-pair scattering. The dominant new effect here 
comes from the renormalization of by the scattering, 
due to U, of a spint^spinj excitation with a Fermi sea 
pair. The corresponding interaction process is described 
by the amplitude $ in Eq. |T]) and is shown schematically 
in Fig. 4(a). The Fermi sea pair (/x, a) is created by in- 
teracting with the spin| electron via U. Such scattering 
gives a contribution oc Uj^av ^aup. to Xf*. I n Fig 3(c) 
we plot this correlation contribution, both for Q close to 
the X-point and for small Q, as function of momentum 
\x for n — 0.6 where the softening is pronounced. We con- 
sider momenta n\\Q (/i?, contribution <&") and momenta 
/i_LQ (fiy, contribution <& ± ). As can be seen in Fig. 3(c), 
the largest correlation contribution comes for /i||Q close 
to the Fermi surface (which for the concentrations of in- 
terest is close to the zone boundary) and for Q close to 
the zone boundary. In Fig3(d) we compare the spinwave 
energy from the full calculation with the results obtained 
by neglecting $- L and/or <P. It is clear that the strong 
softening of the spinwave dispersion as compared to the 
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FIG. 4: (Color online) (a) Schematic describing the scatter- 
ing of spinj-^spinj. electronic excitation with Fermi sea pair 
(/Lt,a) due to U. (b) J c (n) for Jaf = 0.012t, U = 25t. For 
J < J c , the ferromagnetic state \F) is not the ground state. 



RPA comes from $H, i.e. from the renormalization of 
by the scattering of a spin|— >spinj. excitation with a 
Fermi sea pair for momenta /i along T-X. 

With decreasing J/t, the magnon energy for interme- 
diate n turns negative at the X-point while the magnon 
stiffness is still positive. This variational result allows us 
to conclude instability of the ferromagnetic state. On the 
other hand, for small n, the spinwave energy first turns 
negative at the (jr, tt) point (antiferromagnetic correla- 
tions). Finally, for larger n, the spinwave energy turns 
negative at small momenta first, D < 0. By identifying 
the minimum values of J, J c {n), where uq > for all 
momenta, we can definetely conclude, due to the varia- 
tional nature of our calculation, that the ground state is 
not ferromagnetic for J < J c . On the other hand, for 
J > J c , the stability of |F) is not guaranteed. 

J c (n) is shown in Fig^b). By comparing to the RPA, 
it is clear that the pair-magnon correlations lead to a 
very pronounced upward shift of the ferromagnetic phase 
boundary. While for large n the correlation effects dimin- 
ish, and the RPA becomes exact at n—1, for n <0.7 the 
deviations from the RPA exceed 100%. As n decreases 
further, the RPA fails completely and we can conclude 
that it grossly overestimates the stability of the ferro- 
magnetism. Even though additional effects (e.g. phase 
separation @,Q and charge ordering [2(|) will further in- 
crease J c (n) for some n, our variational calculation allows 
us to conclude that Fermi sea pair-magnon correlations 
are strong in the manganites and should be treated be- 
yond the mean field theory of Refs. (R liol]. 

We conclude that non-perturbative long range 
electron-hole pair-magnon correlations play a very im- 
portant role in the spin dynamics of the manganites. 
Most important is the strong softening of the spinwave 
dispersion and the decrease in the stability of the ferro- 
magnetic state. These correlation effects depend sensi- 
tively on the onsite Coulomb repulsion and on its inter- 
play with the magnetic exchange and superexchange in- 
teractions. We propose that the scattering of magnons by 



charge excitations plays an important role in interpret- 
ing recent experiments [11] . Our work can be extended 
to other itinerant ferromagnetic systems (e.g. III(Mn)V 
semiconductors) that are far from the strong coupling 
limit. The correlations discussed here should also play 
an important role in the ultrafast magnetization dynam- 
ics measured by pump-probe optical spectroscopy [21 1 . 
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